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We present theoretical and experimental preparations for an indirect search for new physics (NP) using the 
rare decay Ba — > K*°fj, + fj,~. We design new observables with very small theoretical uncertainties and good 
experimental resolution. 



1. Introduction 

At the start of the LHC we are confronted with 
the experimental fact that all data on flavour ob- 
servables from Babar, Belle, CLEO and also from 
DO and CDF are consistent with the Standard 
Model (SM) predictions [TJ. This implies that 
generic new physics (NP) contributions in K — K 
mixing for example guide us to a new-physics 
scale of 10 3 - 10 4 TeV depending if the new con- 
tributions enter at loop- or tree-level. This is in 
strong contrast to the working hypothesis of the 
LHC that there is NP "around the corner" at 1 
TeV in order to stabilise the Higgs boson mass. 
Therefore, any NP at the 1 TeV scale has to have 
a non-generic flavour structure and we have to un- 
derstand why new flavour-changing neutral cur- 
rents (FCNC) are suppressed. Rare decays and 
CP violating observables allow an analysis of this 
flavour problem. 

The crucial problem in the new physics search 
within flavour physics is the optimal separation of 
NP effects from hadronic uncertainties, ft is well 
known that inclusive decay modes are dominated 
by partonic contributions; non-perturbative cor- 
rections are in general rather small |2l3j . Also 
ratios of exclusive decay modes such as asymme- 
tries are well suited for the new-physics search. 
Here large parts of the hadronic uncertainties 
partially cancel out; for example, there are CP 
asymmetries that are governed by one weak phase 
only; thus the hadronic matrix elements cancel 
out completely. It is the latter opportunity which 
represents the general strategy followed by LHC b 



K*°(i + [i was first 



for the construction of theoretically clean observ- 
ables. 

In this letter we briefly discuss the theoretical 
and experimental preparations for an indirect NP 
search using the rare decay -B^ — > K*° /i + fi~ based 
on the QCDf/SCET approach g]. QCD correc- 
tions arc included at the next-to-leading order 
level and also the impact of the unknown A/mb 
corrections is made explicit 

The exclusive decay -Bd - 
observed at Belle [5] . It offers a rich phenomenol- 
ogy of various kinematic distributions beyond the 
measurement of the branching ratio. We note 
that some experimental analyses of those angu- 
lar distributions are already presented by the 
B factories [6 7 9 10J . Those experimental re- 
sults already have a significant impact on the 
model-independent constraints within the mini- 
mal flavour violation approach [8] . 

Large increase in statistics at LHC b [11|12|13] 
ioi'Bd — > A'*°/i + /i~ will make much higher preci- 
sion measurements possible. There are also great 
opportunities at the future (Super-)_B factories in 
this respect |14ll5ll6ll7j . 

Previously proposed angular distributions and 
CP violating observables in£> d — > K*°[i + [i~ are 
reviewed in Ref. [53], and more recently QCDf 
analyses of such angular distributions [24 25j and 
CP violating observables [26], based on the NLO 
results in Ref. [27], were presented. 
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2. QCD factorization, SCET 

Regarding the hadronic matrix elements of ex- 
clusive modes, the method of QCD-improved 
factorization (QCDf) has been systemized for 
non-leptonic decays in the heavy-quark limit. 
This method allows for a perturbative calculation 
of QCD corrections to naive factorization and is 
the basis for the up-to-date predictions for exclu- 
sive rare B decays in general [TH] . 

A quantum field theoretical framework was 
proposed - known under the name of soft- 
collinear effective field theory (SCET) - which 
allows for a deeper understanding of the QCDf 
approach |19l20j . In contrast to the heavy-quark 
effective theory (HQET), SCET does not corre- 
spond to a local operator expansion. HQET is 
only applicable to B decays, when the energy 
transfer to light hadrons is small, for example to 
B — > D transitions at small recoil to the D meson. 
HQET is not applicable, when some of the out- 
going, light particles have momenta of order m&; 
then one faces a multi scale problem that can be 
tackled within SCET. 

There are three scales: a) A — few x Aqqd the 
soft scale set by the typical energies and momenta 
of the light degrees of freedom in the hadronic 
bound states; b) nib the hard scale set by the 
heavy-6-quark mass and also by the energy of 
the final-state hadron in the -B-meson rest frame; 
and c) the hard-collincar scale /ih c = y/rribA ap- 
pears through interactions between soft and en- 
ergetic modes in the initial and final states. The 
dynamics of hard and hard-collinear modes can 
be described perturbatively in the heavy-quark 
limit mfc — * oo. Thus, SCET describes B decays 
to light hadrons with energies much larger than 
their masses, assuming that their constituents 
have momenta collinear to the hadron momen- 
tum. 

However, we emphasize that within the 
QCDf/SCET approach, a general, quantitative 
method to estimate the important A/mb correc- 
tions to the heavy-quark limit is missing which 
has important phenomenological consequences. 

A careful choice of observables needs to be 
made to take full advantage of the exclusive decay 
Bd — > K*° u, + [i~ , as only in certain ratios such 



as CP and forward-backward asymmetries, the 
hadronic uncertainties cancel out making such ra- 
tios the only observables that are highly sensitive 
to NP. 

Within the QCDf/SCET approach one finds 
crucial form factor relations [21j which simplify 
the theoretical structure of various kinematical 
distributions such that, at least at the leading 
order (LO) level any hadronic uncertainties can- 
cel out. A well-known example of this is the 
zero-crossing of the forward-backward asymme- 
try. In 4] new observables of this kind in the 
Bd — > K*° /i + fi~ decay were proposed which have 
very small theoretical uncertainties and good ex- 
perimental resolution. The only difference to the 
forward-backward asymmetry is that within these 
new observables the hadronic form factors cancel 
out for all values of the dilepton mass. 

3. Theoretical preliminaries 

The decay^ -> E*°£+£~ with K*° -> K~n+ 
on the mass shell is completely described by four 
independent kinematic variables, the lepton-pair 
invariant mass squared, q 2 , and the three angles 
9i, 9k, <f>- Summing over the spins of the fi- 
nal particles, the differential decay distribution 
ofB d -> K* £ + e~ can be written as |28I29I30I31| : 



with 



I = 



h + h cos 29i + I 3 sin 2 9i cos 2cj) 
+I4 sin 29i cos <fi + I5 sin 9i cos <f> 
+Iq cos 9i + If sin 9i sin cj) 
+1$ sin 29i sin <f> + Ig sin 2 9i sin 2(f). 



(1) 



The Ii depend on products of the seven complex 
K* spin amplitudes, A ±L/R , A\\ L / R , A 0L/R , A t 
with each of these a function of q 2 ; the explicit 
formulae are given in the appendix. A t is related 
to the time-like component of the virtual K*, 
which does not contribute in the case of massless 
leptons and can be neglected if the lepton mass 
is small in comparison to the mass of the lepton 
pair. We will consider this case in our present 
analysis. 
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The six complex K* spin amplitudes of the 
massless case are related to the well-known helic- 
ity amplitudes (used for example in 29 30 32] ): 



Ami = (H +lT H. 



2, A = H . 



(2) 



The crucial theoretical input we use in our 
analysis is the observation that in the limit where 
the initial hadron is heavy and the final meson 
has a large energy [21) the hadronic form factors 
can be expanded in the small ratios Aqcd / m>b 
and Aqcb /E, where E is the energy of the light 
meson. Neglecting corrections of order 1 / mt and 
a s , the seven a priori independent B — > K* form 
factors reduce to two universal form factors £j_ 
and £|| [21122] and one finds that the spin ampli- 
tudes at leading order in 1 /m& and a s have a very 
simple form: 



A 



J_L,R 



V2Nm B (l - s) x 



U(Ek*), 



\\L,R 



= -V2Nm B (l 



s) x 



2m 6 (ee) '(cff)x 



(c { 9 cS) tC 10 ) 



Z±(Ek*), 



A 



Nmt 



QL.R 



S\2 



2 m 



= (1-S) 



A'* 



2m 6 (4 eff) -C/ cn; ) 



'(eff)> 



(4° ff) TC 10 ) 
£|| (3) 



with s = q /m B , rhi = wii/ms- Here we ne- 
glected terms of 0(rh 2 K ,). It is important to 
mention that the theoretical simplifications are 
restricted to the kinematic region in which the 
energy of the K* is of the order of the heavy 
quark mass, i.e. q 2 <C m 2 B . Moreover, the influ- 
ences of very light resonances below 1 GeV ques- 
tion the QCD factorization results in that re- 
gion. Thus, we will confine our analysis of all 
observables to the dilepton mass in the range 
lGeV 2 < q 2 < 6GeV 2 . 



4. Construction of theoretically clean ob- 
servables 

By inspection one finds that the distribution 
functions Ii in the differential decay distribution 
(see Eq. (I12P 1 ) are invariant under three symmetry 
transformations which are given explicitly in the 
appendix (see Eqs. ([13M15P ). This implies that 
only 9 of the 12 K* spin amplitudes are indepen- 
dent and that they can be fixed by an full angular 
fit to the 9 independent coefficients of the differ- 
ential decay distribution. Another direct conse- 
quence is that any observable based on the differ- 
ential decay distribution has also to be invariant 
under the same symmetry transformations. 

Besides this mandatory criterium there are fur- 
ther criteria required for an interesting observ- 
able. [Simplicity:] A simple functional depen- 
dence on the 9 independent measurable distribu- 
tion functions; at best it should depend only from 
one or two in the numerator and denominator of 
an asymmetry. [Cleanliness:] At leading order 
in A I nib and in a s the observable should be inde- 
pendent of any form factor, at best for all q 2 . Also 
the influence of symmetry-breaking corrections at 
order a s and at order A/ nib should be minimal. 

[Sensitivity:] The sensitivity to the cff ^ Wil- 
son coefficient representing NP with another chi- 
rality than in the SM should be maximal. [Pre- 
cision:] The experimental precision obtainable 
should be good enough to distinguish different NP 
models. 

In the limit where the K*° meson has a large 
energy, only two independent form factors occur 
in A 0L/R and in A X l/r and A\\ L/R . Clearly, any 
ratio of two of the nine measurable distribution 
functions proportional to the same form factor 
fulfil the criterium of symmetry, simplicity, and 
theoretical cleanliness up to A/ nib and a s correc- 
tions. However, the third criterium, a sensitivity 
to a special kind of NP and the subsequent re- 
quirement of experimental precision, singles out 
particular combinations. In [3] we focused on 
new right-handed currents. Other NP sensitiv- 
ities may single out other observables as will be 
analysed in a forthcoming paper [33) . 
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Figure 1. For A T , theoretical errors (top), ex- 
perimental errors (bottom) as a function of the 
squared dimuon mass, see text for details. 



5. Results 

The first surprising result is that the previously 
proposed quantity , 

A (D r_-r + _ = g(Vl) 

with r± = \H^\ 2 + \H^\ 2 does not fulfil the 
most important criterium of symmetry while it 
has very attractive new physics sensitivity [24l25j . 
Therefore, it is not possible to extract A)p from 
the full angular distribution which is constructed 
after summing over the spins of the final parti- 
cles. Because it seems practically not possible to 
measure the helicity of the final states on a event- 
by-event basis, A)p cannot be measured at either 
LHC& or at a Super- B factory with electrons or 
muons in the final state. 

One finds that the well-known quantities, 
the forward-backward asymmetry Ayr and the 



Figure 2. Aj, , as in Fig.l. 
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Figure 3. A^\ as in Fig.l. 
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longitudinal K* polarization Fl fulfill the sym- 
metry but they include larger theoretical uncer- 
tainties due to the fact that the form factors do 
not cancel at leading order level for all dilepton 
masses. Moreover, the sensitivity to right-handed 
currents is marginal as it is shown below, 

FB 2 |A p + \A {l \2 + \A ± \* {0) 
where for i, j = 0, ||, _L 

MA) = A lL (q 2 )A* L (q 2 ) + A m (q 2 )A; R (q 2 ), 



2 



Fi(<?) IAoP + IAhP + IAJ 2 ' (6) 

In contrast, the following three observables, 
,( 2 ) _ I^l 2 -I^l|l 2 (7] 

A (3) = \ A 0L A 1\L + A 0R A \\r\ (8) 
T ' y/\A \ 2 \A ± \* ' 

^(4) _ \ a olA* ±l - A* R A ±R \ 
T \A* 0L A l[L + A 0R Al R \ ' 

are theoretically clean for all dilepton masses and 
also show a very high sensitivity to right-handed 
currents. 

In the following figures the results on the ob- 
servables, Fl, Afb, A^\ A)p , and A^ are il- 
lustrated: For all the observables the theoretical 
sensitivity is plotted on the top of each figure. 
The thin dark line is the central NLO result for 
the SM and the narrow inner dark (orange) band 
that surrounds it corresponds to the NLO SM un- 
certainties due to both input parameters and per- 
turbative scale dependence. Light grey (green) 
bands are the estimated A/nib ± 5% corrections 
for each spin amplitude while darker grey (green) 
ones are the more conservative A/ mi, ± 10% cor- 
rections. The curves labelled (a)-(d) correspond 
to four different benchmark points in the MSSM 
for righthanded currents (for more details see [4]). 
The experimental sensitivity for a dataset corre- 
sponding to 10 fb _1 of LHC6 data is given in each 
figure on the bottom, assuming the SM. Here the 
solid (red) line shows the median extracted from 
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the fit to the ensemble of data and the dashed 
(black) line shows the theoretical input distribu- 
tion. The inner and outer bands correspond to 
la and 2<T experimental errors. 

The observables AjP and A^p offer sensitivity 
to the longitudinal spin amplitude Aqlr m a con- 
trolled way compared to the old observable Ft'. 
the dependence on both the parallel and perpen- 
dicular soft form factors £||(0) and £j_(0) cancels 
at LO. A residual of this dependence may appear 
at NLO, but as shown in Figs. [2]and[3l it is ba- 
sically negligible. It is also remarkable that for 
Ay and A$p at low q 2 the impact of this uncer- 
tainty is less important than the uncertainties due 
to input parameters and scale dependence. The 
observables A^f 1 and A^ 1 also present a different 
sensitivity to C 7 via their dependence on A l : r 

(2) 

compared with A T . This may allow for a par- 
ticularly interesting cross check of the sensitivity 
to this chirality flipped operator 7 ; for instance, 
new contributions coming from tensor scalars and 
pseudo-scalars will behave differently among the 
set of observables. 

Another remarkable point that becomes clear 

(2) 

when comparing the set of clean observables A T , 
Aa? and A^ versus the old observables Ft and 
Afb concerns the potential discovery of NP, in 
particular of new right-handed currents. There 
are large deviations from the SM curve from 
the ones of the four supersymmetric benchmark 

(2) 

points. A large deviation from the SM for A T \ 

A^ or A^ can thus show the presence of right- 
handed currents in a way that is not possible with 
Ft or Afb- In the latter cases the deviations 
from the SM prediction of the same four repre- 
sentative curves are marginal. 

In the experimental plots we find a good agree- 
ment between the central values extracted from 
the fits and the theoretical input. Any deviations 
seen are small compared to the statistical uncer- 
tainties. The experimental resolution for Ft is 
very good but with the small deviations from the 
SM expected this is not helpful in the discovery of 
new right-handed currents. Comparing the theo- 
retical and experimental figures for the other ob- 
servables it can be seen that in particular A T 
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Figure 6. Belle (black/blue) and BaBar 
(grey/red) data points on Ft and on Afb with 
SM predictions and weighted SM averages over 
the bin q 2 S [1 GeV 2 , 6 GeV 2 ] 



show great promise to distinguish between NP 
models. 

Finally, let us mention that the old observ- 
ables Ft and Afb are already accessible to the 
BaBar [TOI34j and Belle 35] experiments. The first 
measurements are shown in Fig. [B] with the SM 
predictions and the weighted SM averages over 
the bin q 2 e [1 GeV 2 , 6 GeV 2 ]. All the present 
data is compatible with the SM predictions. For 
example, the first measurement of the Babar col- 
laboration on Ft in the low-q 2 region is given as 
an average over the bin q 2 £ [4m 2 , 6.25 GeV 2 ]: 

F L ([4m 2 ,6.25GeV 2 ]) = 0.35 ± 0.16 ± 0.04; (10) 

while the theoretical average, weighted over the 
rate, using the bin, q 2 e [1 GeV 2 , 6 GeV 2 ], based 
on our results is given by: 

F L ([lGeV 2 ,6GeV 2 ]) -0.86^. (11) 
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Here, one should keep in mind that the spectrum 
below 1 GeV 2 is theoretically problematic due to 
the influence of very light resonances; moreover 
the rate and also the polarisation Fl are chang- 
ing dramatically around lGeV 2 . Therefore, we 
strongly recommend to use the standard bin from 
1 GeV 2 to 6 GeV 2 in all future measurements. 

6. Summary 



The full angular analysis of the decay Bd — > 
at the LHC& experiment offers great 
opportunities for the new physics search. New 
observables can be designed to be sensitive to a 
specific kind of NP operator within the model- 
independent analysis using the effective field the- 

(2) (3) 

ory approach. The new observables A y T , A y T ' 
and A^p are shown to be highly sensitive to right 
handed currents. Clearly, theoretical progress on 
the A/rrib corrections would enhance their sensi- 
tivity significantly and would be highly desirable 
in view of a possible upgrade of the LHC& exper- 
iment. Moreover, we have shown that the previ- 
ously discussed angular distribution A^p cannot 
be measured at either LHC & or at a Super- B fac- 
tory. 



h = 



h = 



-|A 0L | 2 cos 2 9 K + (L^R) 
csin 2 9k + dcos 2 9k, 
1 
2 



(|A ±i | 2 -|A|| L | 2 )sin 2 9 K + (L^R) 



= esin 9k, 
h = 1= %l(A 0L A\ L ) sin 29 K + (L - R) 

= f sin 29 K , 
h = V2 $l(A 0L A* ±L ) sin 29 K - (L -» R) 
= g sin 29k, 

^||L^li)sin 2 9 K -(L^R) 



hsin 9k, 



It = V2 



Z(A OL A; L )sm29 K -(L^R) 



j sin 20 K 
1 



<3{A 0L Al L )sm26 K + (L ^ R) 



V2 
k sin 29k, 

3(^ L A ±L )sin 2 9 K + (L^R) 
to sin 2 9k- 



(12) 
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Appendix 

We add here the explicit formula for the distri- 
bution functions and their symmetries: 

In the massless limit, the distribution functions 
Ii depend on products of the six complex K* spin 
amplitudes, A ±L/R , A\\ L/R , A 0L/R : 

h = \{\A^ L \ 2 + \A\ ]L \ 2 + (L^R)) sin 2 9 K + 
+ (|,4oL| 2 + |Am| 2 )cos 2 9 K 
= a sin 2 9k + b cos 2 9k , 

h = i(|^L| 2 + |^|| L | 2 )sin 2 ^ + 



Taking into account a = 3c and b = —d, we are 
left with 9 independent parameters which can be 
fixed experimentally in a full angular fit. 

The distribution functions are invariant un- 
der the following three independent symmetry 
transformations of the spin amplitudes as one 
easily verifies, using the explicit formulae given 
above: (1) a global phase transformation of the 
L-amplitudes 



A' ±L 


= e^Mxi 


A \\L 


= e^A llL , 


AqL 


= e^A 0L] 



(2) a global transformation of the i?-amplitudcs 





= e* 


bR A ±R 


A \\R 


= e" 




Am 


= e u 


iR A 0R ; 
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and (3) a continuous L <-> R rotation 





= + COS0 Axl-\ 


-sm6A* ±R 


A±R 


= -8m6A* XL -\- 


cos 0Aj_r 




= +cos9A 0L - 


sm6A* m 


Aqr 


= + sin 9A* 0L + 


COS OAqr 


A \\L 


= + cos6'v4|| i - 


sin 0A\ R 


A \\R 


= +sin6AT lL + 


cos 0A\\ R . 
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